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When a linear chain of plasmonic nano-particles is subject to longitudinal magnetic field, it 
exhibits optical Faraday rotation. If the magnetized nano-particles are plasmonic ellipsoids arranged 
as a spiral chain, the interplay between the Faraday rotation and the geometrical spiral rotation 
(structural chirality) can strongly enhance non-reciprocity. This interplay forms a waveguide that 
permits one way propagation only, within four disjoint frequency bands; two bands for each direction. 
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Linear chains of identical and equally spaced 
micro-particles have been studied in a number 
of publications It has been shown that 

they allow the propagation of optical modes 
with relatively low attenuation and with no ra- 
diation to the free space. This property is ob- 
tained if the inter-particle distance is smaller 
then the free space wavelength A, and then the 
total width of the modes can be much smaller 
than A. Hence the name "Sub-Diffraction 
Chains" (SDC). SDCs are potential candidates 
for dense integration of optical systems, and 
were proposed asguiding structures, junctions, 
and couplers Very recently, a spiral SDC 

was suggested as a chiral waveguide 

One way waveguides are another type of key 
components in optical systems. Beside their po- 
tential role in reducing disorder effects and un- 
wanted couplings 0, they are used in practice 
as optical isolators and circulators. 

We suggest the SDCs as candidates for one- 
way waveguides. The propagation in plas- 
monic particles SDCs subject to a longitudi- 
nal magnetic field, is non-reciprocal as it pos- 
sesses Faraday rotation p| Q . We suggest an 
SDC in which the plasmonic particles are pro- 
late ellipsoids arranged in a spiral structure as 
in Q, but under longitudinal magnetic field. 
See Fig. [1] We show that the interplay between 
the Faraday rotation and the geometrical spi- 
ral rotation (structural chirality) can strongly 
enhance non-reciprocity, and practically forms 
a one-way waveguide. This one-way waveguide 
is more compact than the one suggested in Q 
(but perhaps more difficult to fabricate), and 
requires an order of magnitude weaker magnetic 
field. At first glance it has a clear physical inter- 



* hadady@eng.tau.ac.il 
' steinber@eng.tau.ac.il 



pretation; a needle-like particle may act locally 
as a polarizer in the sense that it is strongly 
excited by fields aligned with its longest axis, 
and weakly excited by fields normal to this axis. 
Therefore the spiral chain of ellipsoids acts as a 
distributed polarizer that rotates together with 
the Faraday-rotated SDC mode, doing in a dis- 
tributed fashion the same job that a couple 
of polarizers do in a conventional polarization- 
dependent isolator. However, this physical in- 
terpretation is somewhat short of capturing the 
rich physics conveyed by the two-type rotations 
interplay. The one-way property exists even 
with "footballs" with axes ratio of 1:2. Further- 
more, a single spiral chain can support one way 
behavior in four disjoint frequency bands simul- 
taneously, two for each direction. The propaga- 
tion in the "forbidden" direction is practically 
scattered from the chain to the free space; it 
becomes a radiation mode. Principally these 
radiation modes can be blocked by a periodic 
structure as in [7]. It should be emphasized 




FIG. 1. A spiral chain of prolate ellipsoids, subject 
to a longitudinal magnetic field. 

that the effects reported here are not restricted 
to prolate-ellipsoid particles. Thin wire parti- 
cles or any other particles with geometry that 
breaks spherical symmetry would convey essen- 
tially the same physics; the heart of the matter 
lies in the interplay between chirality and Fara- 
day rotation, and not in the specific form of the 
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ellipsoids. The choice of ellipsoids is made to 
simplify the mathematical analysis. Structures 
as in Fig. [T] were fabricated using force medi- 
ating polymer to bend silicon nanopillars [Toj ] . 
that can be coated by metal. Bending angle can 
reach 50 — 60° and it can be tuned by varying, 
e.g., electron beam exposure. A linear array 
of vertically grown silicon nanopillars was de- 
formed by bending the pillars sideways, with 
bending angle that gradually increases along 
the array, thus producing a chain very similar 
to Fig. □ (see Figs. 2,6 in [l(|). 

To study the system, we use the Discrete 
Dipole Approximation (DDA) and polarizabil- 
ity theory. These are standard tools used in 
many works on SDCs They hold when 

the particle diameter D is much smaller than 
the wavelength, and the inter-particle distance 
d is large compare to D. Studies show excel- 
lent agreement with exact solutions even when 
d = 1.5D [ll|. Finally note that the chain is 
periodic only for rational A9/2ir, but in a refer- 
ence frame that rotates together with the spiral 
it is periodic for any AO. Our analysis does not 
assume a' priori any periodicity. 

If a small particle with electric polarizability 
a is subject to an exciting electric field whose 
local value in the absence of the particle is E L , 
its response is described by the electric dipole 
p = aE L . The tensor-polarizability of a gen- 
eral ellipsoid made of an anisotropic material 
e can be found in 12] for the static case. In 
the dynamic case it needs to be augmented to 
incorporate radiation loss [l3j . If the ellipsoid 
principal axes are aligned with the x, y, z axes, 
its dynamic polarizability a is obtained via 
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a x = 2a has N x w 0.1736, N y « 0.4132. A nee- 
dle with a x = 10a has N x rj 0.02, N y « 0.49. 

If the plasmonic particle is subject to an ex- 
ternal DC magnetic field Bo = zBq, e is that 
of a magnetized plasma [l4[. Then 
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where uj v and ajf, = —q e B$/m e are the plasma 
and cyclotron frequencies. Equations (HJ-Q 
describe the ellipsoid at the origin. The n-th 
particle polarizability, a n , is 
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where T„ is a rotation matrix, whose non-zero 
entries are t\\ = £22 = cosnAf?, £33 = 1, £12 — 
—t-ii = sinnAd. The electric field at (0,0,2;) 
due to a short dipole p at (0, 0, z') is given by 
the matrix relation 



E(z) = eo 1 A(z-z')p 
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here Ai = diag(l,l,0), A 2 = diag(-l, -1, 2). 
We express now the local exciting field of the 
m-th particle in the chain as a sum of contri- 
butions from all its neighbor-dipoles, and apply 
a m . The result relates the m-th dipole excita- 
tion p m to its neighbors 
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Here k is the free space wavenumber. I3 is 
the 3x3 identity matrix; together with the 
factor — it represents the particle radia- 
tion loss [13j . a h is the normalized Hermi- 
tian matrix representing the ellipsoid geome- 
try and material. V = Aita x a y a z /?> is the el- 
lipsoid volume and a x ,a y ,a z are its semiaxes. 
L = diag(A x , N y ,N z ) is the depolarization ma- 
trix whose entries are obtained by elliptic inte- 
grals and satisfy ^2 U N u = 1 [12] . To get a feel- 
ing of the numbers involved, note that a sphere 
of radius r has a u — r, N u — 1/3. A prolate 
ellipsoid (a x > a y = a z = a) 
(l-A^/2, A =(l-e 2 )[ln(l±f) 

where e = y/l — a 2 ja 2 x . A 



has N v = N z = 
2e]/(2e 3 ), 
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Note that T_ r „ 
Also, T„ and A commute: T„A(z) = A(z)T r , 
Hence we obtain the shift-invariant difference 
matrix-equation 
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where p r n is the rotated dipole 
Pn = T„p n . 
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For infinite chains, every solution of ([§]) can al- 
ways be expressed as 
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Substituting it back into Eq. ([9]), using 
Eqs. JI])-© and rearranging, we get 

(a^-iI 3 -C)pS=0, (12) 

where C is a matrix defined by the sum 

C = PT E Mnd)T n e~^ nd . (13) 
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A non-trivial solution for the vector exists 
iff the determinant vanishes. In addition, the 
structure supports an optical mode that does 
not radiate to free space iff the radiative-loss 
part represented by — 1L3 cancels out, rendering 
/3 real j4|. This can be achieved only with the 
help of C since 1 has a real diagonal. Hence, 
the guided mode satisfies simultaneously 

Dot (c^ 1 - il 3 - C) = 0, (14) 
-Im(diag(C))=I 3 . (15) 

These equations should be solved for the dis- 
persion curve f3(u>). Where, in the parameter 
space d, A9, j3, u> should we look for the solu- 
tion? Equation © observes the system within 
a reference frame that rotates together with 
the spiral. In this frame the structure is pe- 
riodic, and from (fTTj) (3 is 2ir/d periodic. When 
observed in the "lab" frame, j3 splits to two 
wavenumbers /3i = /3±A9/d. This is easily veri- 
fied: apply T_„ to Eq. (fTTj) . and use T_„ eigen- 
values and eigenvectors Ai,2 = exp(±mA#) and 
v\ t 2 — (1, ±z, 0) / V2 (A3,«3 are ignored as we 
look for transverse excitations). The result is 

Pa = (^*-pS)«ie m( ^ +Ae) + K-pS)^e m( ^- Ae) . 

(16) 

In the lab frame, any SDC guided mode must 
have a real wavenumber within the interval 
k < \/3e\ < n/d Q, establishing our search do- 
main kd < \(3d±A8\ < it. The entries of C 
can be evaluated efficiently using the polyloga- 
rithm functions as in Q, with appropriate ±A9 
shifts (we skip the details as they pertain only 
to numerical efficiency). Finally, note that a 
of a prolate ellipsoid has two resonant frequen- 
cies, due to the two different axes. Hence, an 
ellipsoids SDC possesses two frequency bands 
of guided modes. As we show below each of the 
bands supports two one-way waveguides. 

We solved ([I4|)-([l5]) for w(/3), with the fol- 
lowing parameters. The ellipsoid axes are a x — 
0.25d,a y — a z = 0.5a x and d = A p /30, where 
\ p is the Wp wavelength. The results are shown 
in Fig. [2] When no magnetic field is present and 
all ellipsoids are lined up {AO = 0), they form 
an SDC with transmission bandwidth Aui = 
3.229 x lCP 3 ^,, around the central frequency 
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FIG. 2. The chain dispersion in the lower mode. 
Blue lines are guided modes and black-dashed lines 
are radiation modes, (a) all ellipsoids are lined up 
{A8 = 0) and H = 0. Note the gap at \0\ < k, 
confined between the two radiation modes. Light 
lines are shown in red. (b) A9 — 70° and Ho ~ 0. 
The gap splits symmetrically, (c) The same as (a) 
but with a magnetic field of uib = 0.01cj p . (d) A 
spiral with A8 = 70° , and magnetic field of cut = 
O.Oluip. The gaps of (b) shift to non-symmetrical 
locations. One-way guiding in +z (~z) direction is 
obtained in point A (B). 



cjoi = 0.41702w p - the lower band, and a sec- 
ond SDC with A0J2 = 2.05 x 10~ 3 cj p around the 
central frequency loq 2 = 0.643052w p - the upper 
band. Both have dispersion properties similar 
to the transverse mode in 0] . For these param- 
eters, the DDA works very well [ll|. The lower 
band is shown in Fig. [2^,. The vertical axis is 
Slo = lj — luq 1 , normalized to Au)% . Note the 
gap at |/3 1 d < k. The dispersion lines shown 
as dashed lines at the edges of this gap arc 
practically touching the light-line cone = k 
(shown in red), representing plane- wave modes 
that practically do not interact with the chain 
[H, and can be considered as radiation modes. 
In Fig. [2b a spiral rotation of Ad = 70° is 
added to the chain, yet with B = 0. The 
gap splits in two symmetrically. For clarity, the 
nearly vertical lines of the radiation modes were 
omitted. In Fig. [2J; the ellipsoids are lined up 
(A8 = 0), with a magnetic field of cot, = O.OlWp . 
The dispersion curve is symmetric, but a solu- 
tion for p reveals Faraday rotation. Finally, 
Fig. [2ji shows the dispersion for A9 = 70° 
and Ub = O.Olajp. The magnetic field shifts 
the two gaps of (b), formed by the spiral ro- 
tation, to non-symmetrical locations. At point 
A (a; = uja, Suj = 0.0866Acj) the chain permits 
propagation only in the +z direction (positive 
group velocity). Formally, the horizontal line of 
point A intersects also the light-cone that rep- 
resents propagation in the opposite direction. 
However this line represents a background radi- 
ation mode that is practically not excited as its 
interaction with the chain is very weak. Like- 
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wise, at point B (w = ujb, Soj — 0.0278Aw) the 
waveguide permits propagation only in the — z 
direction. Two one-way guiding are supported, 
one for each direction. 

Essentially the same holds also for the SDC 
modes in the upper band, where we have Slo = 
0.1944Aw for point A, allowing only — z propa- 
gation, and 5u> = O.lllAa; at point B, allowing 
only +z propagation. This is shown in Fig. [3] 
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FIG. 3. The same as Fig.[2]but for the higher mode. 
Note the anti-symmetry with respect to the lower 
mode: one-way —z (+z) is obtained in point A (B). 



cally applied. Hence we use the DDA also here. 
Then Eq. (|9]) still applies, but with two mi- 
nor changes. The summation becomes finite so 
the shift-invariance is lost. Also, the excitation 
of the forced particle at the origin is known: 
Pq = x. Thus Eq. ([9]) becomes a finite ma- 
trix equation for N — 1 unknowns p r n , n ^ 0, 
that is solved numerically. FigureHk shows \p\-\ 
on a logarithmic scale for excitation frequencies 
uja us at the lower band, and Fig. [4]o shows 
the same but at the higher band. In all cases, 
propagation in the "forbidden" direction decays 
within a distance of O(X), becoming at least 
two orders of magnitude weaker. Since the sim- 
ulated chain is of finite length, a reflected mode 
is excited whenever the propagating mode hits 
the chain end. This mode interferes with the 
propagating mode, and generates strong oscil- 
lations. However, since the reflected mode is 




To verify these properties, we have simulated 
a finite chain of N — 800 ellipsoids with the 
above parameters. The ellipsoid at the origin is 
forced by a unit amplitude. When hundreds of 
deep sub-wavelength particles are considered, 
full 3D solutions based e.g. on FDTD tech- 
nique would require tens of thousands of grid- 
points and may become too heavy to handle; 
an equivalent model that encapsulates the main 
physics, with considerably fewer unknowns is 
called for. Recall that for the parameters used 
here the DDA is usually in excellent agree- 
ment with exact full-wave 3D numerical solu- 
tions 11], whenever the latter can be practi- 
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FIG. 4. The chain response to a unit excitation of 
the ellipsoid at the origin, (a) Lower mode. Curves 
A, B correspond to uj — loa,^>b, respectively, (b) 
Upper band. 



propagating in the "forbidden" direction, it de- 
cays exponentially, and so do the oscillations. 
This is seen clearly in the figures. 
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